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Concise theory of chiral lipid membranes 
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A theory of chiral lipid membranes is proposed on the basis of a concise free energy density which 
includes the contributions of the bending and the surface tension of membranes, as well as the chi- 
rality and orientational variation of tilting molecules. This theory is consistent with the previous 
experiments [J.M. Schnur et al., Science 264, 945 (1994); M.S. Spector et al, Langmuir 14, 3493 
(1998); Y. Zhao, et al, Proc. Natl. Acad. Sci. USA 102, 7438 (2005)] on self- assembled chiral 
lipid membranes of DC8,gPC. A torus with the ratio between its two generated radii larger than y/2 
is predicted from the Euler-Lagrange equations. It is found that tubules with helically modulated 
tilting state are not admitted by the Euler-Lagrange equations, and that they are less energetically 
favorable than helical ripples in tubules. The pitch angles of helical ripples are theoretically esti- 
mated to be about 0° and 35°, which are close to the most frequent values 5° and 28° observed in 
the experiment [N. Mahajan et al., Langmuir 22, 1973 (2006)]. Additionally, the present theory can 
explain twisted ribbons of achiral cationic amphiphiles interacting with chiral tartrate counterions. 
The ratio between the width and pitch of twisted ribbons is predicted to be proportional to the rela- 
tive concentration difference of left- and right-handed enantiomers in the low relative concentration 
difference region, which is in good agreement with the experiment [R. Oda et al., Nature (London) 
399, 566 (1999)]. 

PACS numbers: 68.15.+e, 87.10.+e, 61.30.-v 



I. INTRODUCTION 

Since tubules were fabricated successfully from chi- 
ral lipid molecules, chiral lipid structures have attracted 
much experimental attention [l|, 0, H, 0, H, H S B 

El- Most of chiral structures in experiments are 
self-assembled from T)C%$PC which is a typical chiral 
molecule. Schnur et al. have observed that the spheri- 
cal vesicles in solution have very weak circular dichroism 
signal while tubules have strong one Q. Spector et al. 
have investigated the chiral lipid tubules formed from 
various proportions of left- and right-handed DCs.gPC 
molecules and found that they are of similar radii 
which reveals that the radii of chiral tubules do not de- 
pend on the strength of the molecular chirality. Fang's 
group has carefully resolved the molecular tilting order 
in the tubules and concluded that the projected direction 
of the molecules on the tubular surfaces departs 45° from 
the equator of the tubules at the uniform tilting state [H . 
Helical ripples in lipid tubules are also observed by the 
same group with atomic force microscopy § • Their pitch 
angles are found to be concentrated on about 5° and 28° 
Cholesterol is another kind of chiral molecules used 
in the experiments where helical stripes with pitch angles 
11° and 54° are usually observed [7, 8]. Additionally, Oda 
et al. have reported twisted ribbons of achiral cationic 
amphiphiles interacting with chiral tartrate counterions 
S El • It is found that the twisted ribbons can be tuned 
by the introduction of opposite-handed chiral counterions 
in various proportions [9(. From the experimental data 
Q , we see that the ratio between the width and pitch of 
the ribbons is proportional to the relative concentration 
difference of left- and right-handed enantiomers in the 
low relative concentration difference region. Can we in- 
terpret all or at least most of above experimental results 



within a unified theory? 

There are several theoretical discussions on chiral lipid 
membranes (CLMs) in the previous literature, where the 
chiral molecules are assumed to be in a Smectic C* phase 
at which the direction of the molecules is tilted from 
the normal of the membranes at a constant angle. The 
possible free energy of CLMs is discussed by Helfrich 
and Prost from symmetry arguments (Tl1 |. Their theory 
has been further developed and applied in many studies 
[3 El 0, E El EE Hi • Nelson and Powers have inves- 
tigated the thermal fluctuations of CLMs [l3[. Selinger 
et al. have discussed tubules with helically modulated 
tilting state and helical ripples in tubules [lj, llfj, llq . 
Komura and Ou-Yang have given an explanation to the 
high-pitch helical stripes of cholesterol molecules [ljj • 
Due to the complicated form of the free energy used in 
these theories [ill, EH, El El El EI Ej], it is almost im- 
possible to obtain the general Euler-Lagrange equations 
corresponding to the free energy. Thus one cannot deter- 
mine whether a configuration, such as a twisted ribbon, a 
tubule with helically modulated tilting state or a helical 
stripe, is a true equilibrium structure or not. This diffi- 
culty was always ignored in previous discussion of both 
tubules with helically modulated tilting state and helical 
stripes [3, El [Mill- Can we confront this difficulty 
and construct a more concise theory of CLMs consistent 
with the experiments, in which we can unambiguously 
say which configuration is a genuine equilibrium struc- 
ture? 

We will address these questions in this paper, which is 
organized as follows: In Sec. |TT1 we introduce a concise 
free energy density of CLMs which includes the contribu- 
tions of the bending and the surface tension of the mem- 
branes, as well as the chirality and orientational variation 
of the tilting molecules. In Sec. lIIIl we present the Euler- 
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Lagrange equations for CLMs without free edges and use 
them to explain experimental data @, H, 0, Q. We pre- 
dict a torus with the ratio between its two generated 
radii larger than y/2, which has not yet been observed in 
the experiments on self-assembled CLMs. In Sec. IIV1 we 
present the Euler-Lagrange equations and boundary con- 
ditions for CLMs with free edges and use them to discuss 
experimental data [7|, [8|, IS)- Sec. [V] is a brief summary 
and discussion. In the Appendixes, we briefly derive the 
Euler-Lagrange equations for CLMs without free edges 
and the Euler-Lagrange equations as well as boundary 
conditions for CLMs with free edges through the varia- 
tional method |19( developed by one of the present au- 
thors. Several mathematical details are also put in the 
Appendixes. 

II. FREE ENERGY DENSITY 

Following the above theories [H GJ, OS Q EE ES El , 
we adopt a concise form of free energy density for a CLM 
which consists of the following contributions. 

(i) The bending and surface energy per area is taken 
as Helfrich's form 20] 

G H = (k c /2)(2H + c f -kK + \, (1) 

where k c and k are bending rigidities, and A the sur- 
face tension, cq is the spontaneous curvature reflecting 
the asymmetrical factors between two sides of the mem- 
brane. H and K are the mean curvature and Gaussian 
curvature of the membrane, respectively, which can be 
expressed as 2H = -(1/R l + 1/R 2 ), K = l/RiR 2 by the 
two principal curvature radii R± and R 2 . The curvature 
energy in Eq. JT]) is invariant under the coordinate rota- 
tion around the normal of the membrane surface, but will 
change under the inversion of the normal if cq ^ 0. 

(ii) The energy per area originating from the chirality 
of tilting molecules has the form [l2j 

G ch = -hr m , (2) 

where h reflects the strength of the molecular chirality 
which usually determines the handedness of CLMs — the 
CLMs with the opposite handedness will be observed in 
experiments if h changes its sign Q ■ Without losing the 
generality, we only discuss the case of h > in this paper. 
T m is the geodesic torsion along the unit vector m at some 
point. Here m represents the projected direction of the 
lipid molecules in the experiments B, H, EL [H, B, H, B and 
chiral tartrate counterions in the experiment [9|, [l(| on 
the membrane surface, respectively. If we take a right- 
handed orthonormal frame {ei, e 2 , e%} with e3 being the 
normal vector of the membrane as shown in Fig. [TJ m 
can be expressed asm = cos <pe\ + sin <^e 2 , where <j> is the 
angle between m and ei. At this frame, the curvature 
tensor can be expressed as a matrix with element a, b, 
and c shown in Appendix [A"l With the curvature tensor, 
the geodesic torsion along m can be expressed as [l9j 

T m = &(cos 2 4> — sin 2 (j)) + (c — a) cos sin 0, (3) 



which is the same as the chiral term in the previous liter- 
ature, such as the last term of Eq. (2) in Ref. [Tl| , E q (3 ) 
in Ref. [H[, and the third term of Eq. (2.1) in Ref. 15]. 
In particular, for the principal frame, the above equation 
is simplified as 

r m = (1/Ri - l/R 2 ) cos</>sm0. (4) 

We will confine <fr to the region (— 7r/2,7r/2] because of 
the relation r m (0 + n) = r m (</)). Moreover, it is easy to 
see that the geodesic torsion along the mirror image of 
m with respect to ei changes its sign because <j) i— > — <p 
under the reflection with respect to ei. Thus this term 
breaks the inversion symmetry in the tangent plane at 
each point of the membrane, which allows us to distin- 
guish the handedness. Additionally, from Eq. (01 we can 
see that the minimum of Eq. © , 

G™ n = -{h/2)\l/R 1 -l/R 2 \, (5) 

is reached when m departs from the principal direction 
at a angle +7r/4 or — ir/4 for a given shape of the mem- 
brane. Here the sign in front of 7r/4 depends on the sign 
of (l/Ri — 1/R 2 ). The larger difference between i?i and 
i?2 is, the larger absolute value of G™" is reached. In 
other words, the chiral term favors saddle surfaces (for 
example, the twisted ribbons in Sec. IIV C| whose two 
principal curvature radii R\ and R 2 have opposite signs. 




FIG. 1: Right-handed orthonormal frame {ei,e2,es} at any 
point in a surface where e3 is the normal vector of the sur- 
face. (A) Surface without boundary curve. (B) Surface with 
boundary curve where t is the tangent vector of the boundary 
curve, and b, in the tangent plane of the surface, is perpen- 
dicular to t. 

(iii) The energy per area due to the orientational vari- 
ation of m is taken as [2l| 

G ow = (fc / /2)[(Vxm) 2 + (V-m) 2 ], (6) 

where kf is a constant in the dimension of energy. This 
is the simplest term of energy cost due to tilting order 
invariant under the coordinate rotation around the nor- 
mal of the membrane surface. V is the 2-dimensional 
(2D) gradient operator on the membrane surface, and 
the 2D cross product "x" gives a scalar. By defining a 
spin connection field S satisfying V x S = K [22|, one 
can derive 

(V x m) 2 + (V ■ m) 2 = (V</> - S) 2 (7) 
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through simple calculations. 

The total free energy density adopted in the present 
paper, G = Gh + G c h + G ov , has the following concise 
form: 

G=^(2H + c ) 2 ~~kK + X-hT m + ^ 2 , (8) 

where v = — S. This special form might arguably 
be the most natural and concise construction including 
the bending, chirality and tilting order, for the given 
vector field m and normal vector field 63. Of course, 
using these two fields, one can also construct more gen- 
eral free energy densities which contain much more terms 
[U H H Q [H US G3- Compared with the general 
form, our special form is simplified in two aspects: (i) 
The contributions of the bending of membranes and the 
orientational elasticity of m are taken as the isotropic 
forms; (ii) There is no additional coupling between m 
and the curvature of membranes except the chiral term 
@, which makes it a concise theory. Although our free 
energy density is not new relative to the general form ap- 
pearing in the previous literature, we derive for the first 
time the corresponding Euler-Lagrange equations with- 
out any assumption on the shapes of membranes before 
doing the variation, and then check which configuration 
solves the Euler-Lagrange equations. More importantly, 
in the following, we will see that most of the experimen- 
tal results can be explained in terms of such a concise 
free energy density. 



III. CHIRAL LIPID MEMBRANES WITHOUT 
FREE EDGES 

CLMs without free edges usually correspond to closed 
vesicles. Here a long enough tubule, where the end ef- 
fect is neglected, is also regarded as a CLM without free 
edges. In this section, we show the Euler-Lagrange equa- 
tions of CLMs without free edges (derived briefly in Ap- 
pendix [B]), and then discuss CLMs in spherical, tubular 
and torus shapes. As we will see, our theoretical results 
are in good agreement with the experiments [2, [H, 0, Q • 



as 

2h(n m - H) - k f V 2 4> = 0, (10) 

and 

2k c V 2 H + k c (2H + c )(2H 2 - c H - 2K) - 2XH + P 
+/i[V • (mV x m) + V x (mV • m)] 
+/c / [( Kv -i?)v 2 -Vv:Ve 3 ] = 0, (11) 

where n m and n v are the normal curvature along the di- 
rections of m and v, respectively. They can be expressed 
as Eqs. (|A6|) and (|A7|1 with the curvature tensor. Phys- 
ically, Eqs. l[TU]) and (fTTj) express the moment and force 
balances along the normal e3 at each point of the mem- 
brane. If h and kf vanish, the above two equations de- 
generate into the shape equation of achir al lip id vesicles, 
which has been fully discussed in Refs. [SImIIIIMIII]- 
Two remarks are necessary concerning Eq. (fTUj) : (i) We 
have selected the proper gauge such that V ■ S = 0, or 
else V 2 should be replaced with V 2 </> — V • S. (ii) For 
closed vesicles different from toroidal topology, the tan- 
gent vector field m will have singular points. In this case, 
the right-hand term should be replaced with the sum of 
(5-function, ^\ <7^(r — r^), where r and represent any 
point and singular point in the CLM, while <ii represents 
the strength of the source or vortex at the singular point 

B. Spherical vesicles 

For spherical vesicles of chiral lipid molecules, r m is 
always vanishing because R\ = i?2- Thus the free en- 
ergy © is independent of the molecular chirality and 
permits the same probability of left- and right-handed 
spherical vesicles existing in solution. Naturally, no evi- 
dent circular dichroism signal would be observed, which is 
consistent with the experiment Of course, we cannot 
exclude the other possible explanation of this experiment 
that the lipid molecules are non-tilting in the spherical 
vesicles Hal. 



A. Euler-Lagrange equations 

The free energy for a closed CLM can be expressed as 

F = jGdA + P J dV, (9) 

where dA is the area element of the membrane and dV 
the volume element enclosed by the vesicle. P is either 
the pressure difference between the outer and inner sides 
of the vesicle or used to implement a volume constraint. 
For tubular configuration, we usually take P = 0. 

Using the variational method [l^| , we obtain the Euler- 
Lagrange equations corresponding to the free energy ^ 



C. Tubules with uniform tilting state 

A tubule is regarded as an cylinder with radius p and 
infinite length. Then H = —l/2p, K m — — cos 2 cj>/p, 
and T m = sin 2<j)/2p, where <fr is the angle between m and 
the equator of the cylinder. <fi is a constant for the tubule 
with the uniform tilting state as shown in Fig. HJa) . Thus 
Eq. (fTO)) requires 

(f> = ±7r/4. (12) 

We should keep <f> = tt/4 because it corresponds to the lo- 
cal minimum of —h J T m dA, which is in good agreement 
with the experiment by Fang's group Q where the pro- 
jected direction of the molecules on the tubular surfaces 
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indeed departing 45° from the equator of the tubules with 
the uniform tilting state. 

For the cylindrical shape, we set P = 0. It is easy to 
find V • (mV x m) + V x (mV • m) and (k v — H)v 2 — 
Vv: Ve3 both to be vanishing for <f> — tt/4. Thus Eq. (fTTj) 
is transformed into 



k c (c 2 p 2 - 1) + 2Xp z = 0. 



(13) 



This equation indicates that the radius p of chiral lipid 
tubules is independent of the strength of molecular chi- 
rality h, which is consistent with the experiment [3j where 
the tubular radius is insensitive to the proportion of 
left- and right-handed DCs,9-PC in the tubule. How- 
ever, we cannot exclude the other possible explanation 
that DC 8 ,g-PC molecules with different handedness are 
not mixed in the tubule such that each tubule in solution 
contains merely the one kind of molecules with the same 
handedness ether left- or right-handed [jj, [l|| . 



(b) 






FIG. 2: Tubules: (a) Uniform tilting state; (b) Helically mod- 
ulated tilting state. Arrows represent the projected directions 
{m} of the tilting molecules on the tubules. 



D. Tubules with helically modulated tilting state 

The concept of tubules with helically modulated tilting 
state was proposed by Selinger et al. [15]. The orienta- 
tional variation of the tilting molecules is assumed to be 
linear and confined to a small range. Here we will directly 
investigate the orientational variation by using the Euler- 
Lagrange equation (|10[) without this assumption. Addi- 
tionally, as mentioned in the Introduction, whether the 
tubules with helically modulated titling state are equilib- 
rium configurations or not is not addressed in Ref. [l5|. 
However, our theory will unambiguously reveal that they 
are not equilibrium configurations. Let s and z denote 
the arc length parameters along the circumferential and 
axial directions, respectively. If cj>, the angle between 
m and the equator of the cylinder, is not a constant, 
Eqs. (|10p and (fTTj) are then transformed into 



and 

+A//H 
= 0. 



> sz ) sin 2(f) - 
l/p 2 )/2p- 



(<Pss - 0zz + Hz^s) COS 2(f)] 

kf[(<t? z -tf)/2p+<t>„/p] 

(15) 



where the subscripts s and z represent the partial deriva- 
tives respect to s and z, respectively. The derivations of 
the above two equations are shown in Appendix ID] 

At the helically modulated tilting state, 4> is invariant 
along the direction of a fictitious helix enwinding around 
the tubule as shown in Fig. [D[b). Let -0 be the pitch 
angle of that helix and apply a coordinate transformation 
(s, z) — ► (£, 77) via C, = s cos ip + z sin?/', T) = — s sin-0 + 
z cos tp, where £ is the coordinate along the helix and 77 is 
the coordinate orthogonal to £. In the new coordinates, 
<t> depends only on r\. Changing variable Q = (f> — tp and 
introducing the dimensionless parameters X = v/p an< i 
h = hp/kf, we transform Eqs. (fTOjl and (fTTj) . respectively, 
into 



and 



Q xx = -h cos 2(6 + ip), 



{k c /k f )(l-clp 2 )-2Xp 2 /k f = R, 



(16) 



(17) 



where R = Q 2 , cos 2ip - 9 XX sin 2ip + 2h{2Q 2 x sin 26 - 
XX cos 26). 

The first integral of Eq. (fTtjl) is 

O 2 = p 2 -hsm2{Q + i)), (18) 



with an unknown constant p. Substituting Eqs. (Tl7; 
and l|18p into Eq. (jTTJ) , we obtain the right-hand side of 



k 



f\ 



(f> zz ) + {h/p) cos 20 = 0, 



(14) 



R = p 2 cos2tp + h(Ap 2 - 1) sin 26 

+ 2/i 2 [cos(4e + 2V>) -sin26sin2(6 + ip)]. (19) 

The necessary condition for validity of Eq. (jTTJ) is that the 
right-hand side R is constant, which, in terms of Eq. (|19p . 
holds if and only if h = for varying 6. Therefore, 
tubules with helically modulated tilting state are not ad- 
mitted by the Euler-Lagrange equation pTjl if h ^ 0. 
In other words, the orientational variation of the tilting 
lipid molecules breaks the force balance along the normal 
direction of the tubule at this state, which might rather 
induce helical ripples in tubules. 



E. Helical ripples in tubules 

Assume now that a tubule with radius p undergoes 
small out-of-planc deformations and reaches a new con- 
figuration expressed as a vector {^(1 + y) cos(s/p), p(l + 
y) sin(s/p), z}, where \y\ <C 1 is a function of s and z. For 
simplicity, we take cq = 0, A = 0, kf ~ k c , and h <C 1 in 
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this subsection. As shown in Appendix [El Eqs. (fT0|) and 
(fTTj) can be transformed into 

p 2 {4> ss +4> zz +4'zVz-4'sys~^y4'ss-yzs) = -hcos2(p (20) 

and 

(fc c /fc /P 2 )[(l + p 2 cU + p 2 9 Z2 ) 2 y - 1/2] 
+ (fc c /fc /P 2 )[(l + p 2 d ss + p 2 d zz )y/2 - 2p 2 y zz ] 
+h[2 (4> 2 Z -4> 2 S + <t> sz ) sin 2<f) + {(j) ss - cj) zz + ty, S ) cos 2(f>] 
+(0 2 - 2 )(y + p 2 y ss - p 2 y 2Z - l)/2 

+0 S 2 + y^s - 2 y^«2 - + ^p 2 y s A a <\> z = o, (21) 

where ^ is the angle between m and the equator of the 
tubule. 




FIG. 3: Small amplitude ripples in a tubule with radius p. 
Arrows represent the projected directions {m} of the tilting 
molecules on the ripples' surface. 

Now we consider helical ripples in the tubule where 
<p and y are invariant along the direction of a fictitious 
helix cnwinding around the tubule as shown in Fig.[3J We 
adopt the same coordinate transformation as the above 
subsection, and let § = <p — (Q+ip), where ip is the pitch 
angle of the fictitious helix and O is governed by Eq. (fT!)]) . 
Then the above equations (|20"f and (|2"Tj) are reduced to a 
matrix equation 

= $, (22) 

with * = {§ lV Y and $ = {0,(V2)(1 - 4/r 2 ) sin26}', 
where the superscript 't' represents the transpose. The 
differential operator £ has four matrix elements: 

C n = ~d 2 /d x 2 , 

C12 = —pcos2ipd/dx — sinipcosipd 2 /dx 2 , 
£21 = pcos2ipd/dx — sin V> cos ipd 2 /dx 2 , 
£22 = p 2 cos 2 ip + p 2 cos 2ipd 4 /dx 4 

+ [2p 2 cos 2ip sin 2 ip + {p 2 - l)cos 2 ^}d 2 /dx 2 , 

where p satisfies 

p 2 cos2ip = k c /kf. (23) 



It is not hard to find a special solution of Eq. (|22]) as 

- _ Ml-4/i 2 ) /cos2(e + V;)\ ( 
*~ 4r ^2sin29 J [ > 

with T = /i 2 (3 cos 2 ip + cos 2^ — 8p 2 cos 2ip sin 2 ip — 
4 / u 2 cos 2 ip + 16/^ 4 cos2V'). Considering Eq. (|23p. we can 
prove r > for k c /kf > 0.133 and hence T > for 
k c — kf. 

Note that Eq. (|2"2"|) contains only the first order terms of 
•& and y, which can also be obtained from the variation of 
the free energy F expanded up to the second order terms 
of $ and y. The dimensionless mean energy difference 
between the tubule with helical ripples and that with 
helically modulated tilting state is expressed as 

r W/p /-, \ 

AF = J f d X , (25) 

where W is the period of the ripples along 77 direction as 
show in Fig. [3l Substituting the solution (|24|) into the 
above equation, we have 

h 2 (l - Au 2 ) 2 f w/p 
AF = - (i gr 4/i J J q sin 2 2Qd X < 0, (26) 

which reveals that a tubule with ripples is energetically 
more favorable than that in a helically modulated titling 
state. Additionally, we can easily prove from Eq. (|23|) 
that the pitch angle obeys ip < 45°. All ripples in tubules 
resolved in the recent experiment Q with atomic force 
microscope have indeed the pitch angles smaller than 45°. 
The tubules with 'helically modulated tilting state' ob- 
served in the experiment [281 ] might also be the tubules 
with helical ripples whose amplitudes are below the ex- 
perimental resolution. In this experiment, all pitch angles 
are also smaller than 45°. 

Considering h <§; 1, we have 6 X ss p from Eq. (fT5|) and 
then 

pW/p w 2tt, (27) 

due to W being the period of the ripples along 77 direction. 
If the pitch angle -0 = 0, the fictitious helix is a circle, 
and Eqs. (|23[) and (|27[) requires 2np/W w y/k~/kj. If 
ip =/= 0, the fictitious helix is indeed a helix satisfying 
W = 2irpsmip. It is then easy to derive 

cos 2^/ sin 2 V ~ k c /kf (28) 

from Eqs. (|2"3")) and (|2~T|) . Therefore, two kinds of ripples 
are permitted in our theory: one has pitch angle about 
0° ; another satisfies Eq. ([28]) , which gives ip ~ 35° for 
k c ~ kf. Our theoretical results are thus close to the 
most frequent pitch angles (about 5° and 28°) of the 
ripples in tubules observed by Fang's group [5J. 

F. Torus, a prediction 

A torus is a revolution surface generated by a cycle 
with radius p rotating around an axis in the same plane 
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(a) 





kf/kc 



FIG. 4: (a) Torus generated by a cycle rotating around an axis 
in the same plane of the cycle, (b) Ratio of the generated radii 
(v), Eq. (I31|l . as a function of k c /kf. 



of the cycle as shown in Fig. Ufa). The revolution ra- 
dius r should be larger than p. A point in the torus 
can be expressed as a vector {(r + p cosy) cos 9, (r + 
pcos (/?) sin#, p sin ip}. As shown in Appendix [Fl Eq. JTI 
is transformed into 



v + cos ip 



[{v + cos (p)(j>, 



<piv> 



vhp 



cos 20 = 0, (29) 



where cf> is the angle between m and the latitude of the 
torus, while 

v = r/p 

is the ratio between two generated radii of the torus. 

It is found that the uniform tilting state <\> — —it /A sat- 
isfies the above equation (f2"9")l and that makes — J hr m dA 
to take the minimum by considering Eq. (|F12[) . As shown 
in Appendix [Fl with — — it/ 4, Eq. (jlip is transformed 
into 

(2k c -k f )/v 2 + k c (c 2 p 2 - 1) + 2{Pp + \)p 2 

4k c clp 2 - Ak c c p - 2hp + 8Xp 2 + 6Pp 3 

H ; cos if 

v 

bkccZp 2 - 8k c c p - 4hp + 10Ap 2 + 3fc f + 6Pp 3 2 
H t, cos if 

+ 2k c clp 2 - 4k c c a p - 2hp + AXp 2 + 2k f + 2Pp 3 3 



= 0. 



(30) 



If v — > oo and P = 0, the above equation degenerates 
into Eq. (|13[) obeyed by a tubule with radius p at uniform 



tilting state. If v is finite, then Eq. (|30|) holds if and only 
if the coefficients of {1, cos if, cos 2 (f, cos 3 f} vanish. It 
follows that 2Ap 2 = fc c (4pc - p 2 c\) - 3k f + 2hp, Pp 3 = 
2k f — 2k c pco — hp and 



kf/k c 



1 - kf/k c 



(31) 



The relation between v, the ratio of the generated radii 
r and p, and kf/k c is sketched in Fig. BJb). T ne ratio v 
increases with kf/k c . Especially, v = \[2 for kf/k c = 0, 
which leads to the Willmore torus of non-tilting lipid 
molecules [2i| H3| • Since this kind of torus was observed 
in the experiments [U HH, [33| , tori with v > \[2 for < 
kf/k c < 1 might also be observed in some experiments 
on CLMs. 



IV. CHIRAL LIPID MEMBRANES WITH FREE 
EDGES 

In this section, we show the Euler-Lagrange equations 
and boundary conditions of CLMs with free edges, and 
then discuss helical stripes and twisted ribbons. As we 
will see, our theoretical result on twisted ribbons is con- 
sistent with the experiment 



A. Euler-Lagrange equations and boundary 
conditions 

Consider a CLM with an free edge as shown in 
Fig. [T](b). Its free energy can be expressed as 



F = / GdA 



7 



ds, 



(32) 



where dA is the area element of the membrane and ds 
the arc length element of the edge. 7 represents the line 
tension of the edge. 

Using the variational method [l9|, as shown in Ap- 
pendix [C] we obtain two Euler-Lagrange equations corre- 
sponding to the free energy ([3"2"]l as 



2h(K r 



H) 



k f V 2 



= 0, 



(33) 



and 



2k c V 2 H + k c (2H + c )(2H 2 - c Q H - 2K) - 2XH 
+/i[V • (mV x m) + V x (mV • m)] 
+k f [{ Kv -H)v 2 - Vv:Ve 3 ] = 0. (34) 

Additionally, the boundary conditions obeyed by the free 
edge are derived as: 



kfV b = 0, 

G + JKg — 0, 

k c (2H + cq) — kn r< 



(ft./2)sin20 = 0, 



(35) 
(36) 
(37) 



kf g — 2k c Hb — h(vt + 4>) sin 20 + kfK n v t = 0,(38) 
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where K n , r g and K g are the normal curvature, geodesic 
torsion, and geodesic curvature of the boundary curve 
(i.e., the edge), respectively, b is in the tangent plane of 
the membrane surface and perpendicular to t, the tan- 
gent vector of the boundary curve, as shown in Fig.[ljb). 
Vb and vt are components of v in the direction of b and 
t, respectively. H\> is the direction derivative of H with 
respect to b. A dot represents the derivative with respect 
to s. is the angle between m and t at the boundary 
curve. 

It should be noted that Eqs. (|33|) and (|34|) are equiv- 
alent to Eqs. (HI and (JHJ with P = 0. Eqs. (|35 |) -|38] l 
describe the force and moment balance relations in the 
edge. Thus they can also be used for a CLM with several 
edges. If h and kf vanish, the above equations (|33[) - 
(f3"5| degenerate into the Euler-Lagrange equations and 
boundary conditions of open achiral lipid bilayers, which 
were fully discussed in Refs. [34l l35l. 13a. [37| . 

B. Helical stripes 

A helical stripe with pitch T and radius p is shown 
in Fig. In terms of the discussion in Sec. IIII D[ we 
can immediately deduce that a helical stripe with modu- 
lated tilting state is not permitted by the Euler-Lagrange 
equations (|33|) and (f34|) because it can be thought of as a 
ribbon wound around a fictitious cylinder. Thus here we 
only have to discuss helical stripes with a uniform tilting 
state, for which we easily obtain 

= tt/4, (39) 
k c (c 2 oP 2 - 1) + 2Xp 2 = 0, (40) 

from the Euler-Lagrange equations. Here is the angle 
between m and the equator of the fictitious cylinder. 




FIG. 5: Helical stripe. Arrows represent the projected direc- 
tions {m} of the tilting molecules on the stripe's surface. 

In the free edge of the helical stripe, we have = v t = 
0, H = -l/2p, = tt/4 - V, Tm = -h/2p, K g = 0, K n = 



— cos 2 tp/ p, T g = sm2t(j/2p, where if) = T/2irp is the pitch 
angle of the helical stripe. Thus the boundary condition 
(|35| is trivial, and Eqs. ([36|) -(|38 ]) can be transformed into 

(k c /2){c -l/p) 2 + X-h/2p = 0, (41) 
k c {c - l/p) + fccos 2 ip/p- (/i/2)cos2i/> = 0, (42) 
7COS 2 V/P = 0. (43) 

If 7 7^ 0, there exists only a trivial solution ip = ir/2 
to the above equations (|40p (14"3")) . Thus our theory does 
not permit genuine helical stripes with free edges with a 
uniform tilting state. 

C. Twisted ribbons 

A twisted ribbon with pitch T and width W is 
shown in Fig. [6l which can be expressed as a a vec- 
tor {u cos <p,u sin ip, atp} with |w| < W/2, \ip\ < oo and 
|a| = T/2n. As shown in Appendix iGl Eq. (f3"3"]) is trans- 
formed into 




FIG. 6: Twisted ribbons: (a) m is perpendicular to the edges; 
(b) m parallels the edges. Arrows represent the projected 
directions {m} of the tilting molecules on the ribbons' surface. 

If we only consider the uniform tilting state, the above 
equation requires <fi = or 7r/2. It is easy to see that 
= minimizes — h J T m dA for a < while = tt/2 
minimizes — h J T m dA for a > from Eq. (|G12[) .Thus 
we should take = for a < and = 7r/2 for 
a > 0. The former case corresponds to Fig. [6]^ a) where 
m is perpendicular to the edges; the latter corresponds 
to Fig. [6{b) where m is parallel to the edges. As shown 
in Appendix [Ul both for = and 7r/2, Eq. ((3"4")) is trans- 
formed into 

k c c a 2 /{u 2 + a 2 ) 2 = 0, (45) 
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which requires cq — for non-vanishing a. Among the 
boundary conditions ([3"5| - ([3"8)) , only Eq. ([3"6"| is nontriv- 
ial, which gives 

A(l + x 2 )a 2 -(h- jx)\a\ + = 0, (46) 

with x = W/2\a\. 

Guided by the experimental data @, we may assume 
h to be proportional to Rd, the relative concentration 



Minimizing it with respect to \a\ and using Eq. (|47| . we 
obtain 

x = /3Rd + 0(R 3 d ), (49) 

with /3 = 3/io/47- This relation reveals that the ratio be- 
tween the width and pitch of the ribbons is proportional 
to the relative concentration difference of left- and right- 
handed enantiomers in the low relative concentration dif- 
ference region. Eq. (|49p fits well with the experimental 
data with parameter (3 = 0.37 as shown in Fig. \7\ 



0.15 



00 




0.0 0.1 0.2 0.3 0.4 0.5 



FIG. 7: Relation between x = W/2\a\ and R d . The solid 
line is the fitting curve x = 0.37Rd and the dots are the 
experimental data in Ref. [9(. 



V. CONCLUSION AND DISCUSSION 

In this paper, we have focused on a concise free en- 
ergy from which the Euler-Lagrange equations of CLMs 
without free edges and Euler-Lagrange equations as well 
as boundary conditions of CLMs with free edges can be 
derived analytically. Their implications for extant and 
future experiments can be summarized as follows. 



difference of the left- and right-handed enantiomers in 
the experiment, i.e., h — h Rd with a constant h . In 
terms of the experimental data, \a\ — > oo for Rd — > 0. 
Thus Eq. (|46[) requires A = and then 

H = (2k + k f x 2 )/2(h R d - 7«)(1 + x 2 ). (47) 

To determine the relation between x and Rd, we minimize 
the average energy per area with respect to \a\ for given 
W. The average energy per area is calculated as 



(48) 

I 

(i) Our theory predicts the same probability for left- 
and right-handed spherical vesicles existing in solution. 
Thus no evident circular dichroism signal should be ob- 
served, which is consistent with the experiment Q. 

(ii) The radius of a tubule at uniform tilting state sat- 
isfies Eq. (|13p . It reveals that the radius is independent of 
molecular chirality, which is consistent with the experi- 
ment @ where the tubular radius is insensitive to the pro- 
portion of left- and right-handed DC8.9PC in the tubule. 
We find that the projected direction of the molecules on 
the tubular surface departs 45° from the equator of the 
tubule with the uniform tilting state, which is in good 
agreement with the experiment Q. 

(iii) Tubules with helically modulated tilting state are 
not admitted by the Euler-Lagrange equation for 
non-vanishing h. The orientational variation of the tilt- 
ing lipid molecules breaks the force balance along the nor- 
mal direction of the tubule. Thus a tubule with helically 
modulated tilting state is not an equilibrium structure 
within our theoretical framework. 

(iv) Helical ripples in tubules are equilibrium struc- 
tures. The pitch angle of helical ripples is estimated as 
about 0° and 35° for kf ~ k c , which are close to the ex- 
perimental values 5° and 28° observed by Fang's group 

(v) Tori with a ratio of generated radii larger than \[2 
are predicted by our theory, Eq. (|3ip , which have not yet 
been observed in the experiments. 

(vi) Helical stripes with free edges at either uniform 
tilting state or helical modulated tilting state are not 
possible equilibrium configurations. 

(vii) Twisted ribbons satisfy Euler-Lagrange equations 
([55]). ([MP and boundary conditions ([55 [l -([55 |) , The ratio 
between the width and pitch of the ribbons is propor- 
tional to the relative concentration difference of left- and 
right-handed enantiomers in the low relative concentra- 
tion difference region, which is in good agreement with 
the experiment @. 

Finally, we have to list a few open problems which 



I 

- (2k - k f )x + (k f - 2h R d \a\)v / l + s 2 arcsinha: + 27|a|(l + x 2 ) 

t = — — 

a 2 y/l 4- x 2 (arcsinhx + xyl + x 2 ) 
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should be addressed in the future work. 

(i) Within the present theory, we cannot give a simple 
explanation for the experiments on CLMs of pure choles- 
terol @,Q where helical stripes with pitch angles 11° and 
54° are usually observed. It has recently been found that 
the cholesterol helical stripes have good crystal structure 
[38| , which is out of the range of our theory. There might 
be two possible ramifications based on the present the- 
ory: One would be to include anisotropic bending effects 
[39j ] in the free energy density ([5]); Another one would 
be to consider a line tension 7 depending on the angle 
between the directions of the tilting and the free edges 
[H in the free energy <[32]>. 

(ii) We cannot yet interpret the twisted ribbon-to- 
tubule transition with increasing the relative concentra- 
tion difference of the left- and right-handed enantiomers 
reported in the recent experiment by Oda's group [ill ]. 
A possible reason is that the parameters except h in our 
theory are independent of the relative concentration dif- 
ference [Hj]. Additionally, the ribbons and tubules ob- 
served in this experiment are usually multi-bilayer struc- 
tures 42]. Then a decoupling effect [HI, [44[ between 
neighbor bilayers might occur. An extended theory in- 
cluding these two factors might be required to investigate 
the mechanism responsible for the transition. 

(iii) We have ignored the effect of singular points in 
CLMs. Although the features of singular points in 2D 
planar films or achiral nematic spherical, torus vesicles 
and other manifolds have been fully investigated in pre- 
vious literature [4fJ, |4g, |4g, |47J, |48| , it remains a chal leng e 
to study the properties of singular points in CLMs [49| ■ 
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APPENDIX A: BRIEF INTRODUCTION TO 
MOVING FRAME METHOD AND EXTERIOR 
DIFFERENTIAL FORMS 

If we take a frame {ei, e2, 63} at any point r on a sur- 
face as shown in Fig. [H then the infinitesimal tangential 
vector at r is expressed as 

dr = LJiei + o7 2 e 2 (Al) 

and the difference of frame between at points r + dr and 
r is denoted as 

de, = uJijej, (i = 1, 2, 3), (A2) 

where w\, LO2, and u>ij = — Wji, = 1,2,3) are I-forms 
[T^ . |5C| . The repeated subscripts in Eq. (|A2[) and the fol- 
lowing contents represent the Einstein summation con- 
vention. With these I-forms, the structure equations of 



a surface can be expressed as [T^, [5(| 

dwi = L0\2 A LU 2 , 

di02 = W21 A wi, (A3) 
doJij = uj ik A uj k j (i, j = 1, 2, 3), 

and 

(E)-(i (A4 » 

where the symbol 'A' expresses the wedge product be- 
tween differential forms and 'gT is the exterior differential 

operator (l9l . |50| . The matrix ^ ^ ^ ^ is called the cur- 
vature matrix which is related to the mean and gaussian 
curvature by 

2H = a + c, K = ac-b 2 . (A5) 

For a unit vector m = cos(/>ei + sin</>e2, the normal 
curvature along the direction of m can be expressed as 

m 

n m = a cos 2 4> + 2b cos 4> sin <p + c sin 2 (f>. (A6) 

The normal curvature along the direction of an arbitrary 
vector v can be expressed as 

k v = (avf + 2bv 1 v 2 + cvj)/v 2 , (A7) 

where v\ and V2 are the components of v in the directions 
of ei and e2- 

In our following derivations, several relations between 
vector forms and differential forms are used frequently. 
For convenience, we list them below. 



(V x m)dA = d(m ■ dr), 


(A8) 


(V • m)dA = d(*m ■ dr), 


(A9) 


T m dA — m ■ de% A m • dr, 


(A10) 


K m dA = — m ■ de^ A *m • dr, 


(All) 


V4> ■ dr = d<f), 


(A12) 


(y 2 <j))dA = d*d<t>, 


(AI3) 


S • dr = — 


(A14) 


(V • S)dA = -d*Lu X2 , 


(A15) 


(V x S)dA = -duj 12 = KdA, 


(AI6) 


v • dr = d(p + UJ12, 


(A17) 


w 2 dA = v • dr A *v • dr, 


(A18) 


(Vv:Ve 3 )dA = d(\ ■ de 3 ), 


(A19) 



where dA — wi A W2 is the area element, m is a unit 
vector and <fi is the angle between m and ei. S is the spin 
connection and v = V(f)— S. * is the Hodge star operator 
p^ . l5l| which satisfies *loi — L02 and *uj2 — — u>\ . 

Using Eqs. (fM]), (fAT]), CM]), (ESJ) and (|AT7|) . we can 
prove 

(k v - H)v 2 = {v{ - v 2 ){a - c)/2+ 2b Vl v 2 , (A20) 
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and 

[V • (mV x m) + V x (mV ■ m)]dA (A21) 
= d[(v2 cos 2(f) — vi sin 2<f))(Ji + (i>i cos 2(f) + v 2 sin 20)^2] , 

where v\ and v 2 are the components of v in the directions 
of ei and e2 . 

We suggest the reader refer first to Ref. [li| and some 
textbook on the calculus with the moving frame method 
before going on if he is interested in the mathematical 
details. In writing the following contents we assume that 
the reader has been familiar with the skill in Ref. [l9[ 
and good at calculus with the moving frame method. 



APPENDIX B: DERIVATION OF THE 
EULER-LAGRANGE EQUATIONS OF CLMS 
WITHOUT FREE EDGES 

Here the CLMS without free edges will be derived 
through the variational method developed in Ref. [l|| 
with the aid of the moving frame method and exterior 
differential forms, which can highly simplify the calculus 
of variation. 



1. Variation respect to <j> 

Assume 8(f> — S. Through simple calculations, we ar- 
rive at 



functional 



and 



2(H-K m )E, 



5{-v 2 dA) = 2dZ A *(d(f) + w 12 ). 



(Bl) 



(B2) 



Combining with the above two equations and using the 
integral by parts and Stokes' theorem, we derive 

5F = J[2h(n m -H)-k f V 2 (f)}EdA+k f <j> S*v-dr, (B3) 

where the second term is the integral along the bound- 
ary curve, which vanishes for CLMS without free edges. 
Since S is an arbitrary function, from SF = we derive 
Euler-Lagrange equation (JTUJ) . We have used the assump- 
tion V ■ S = when we write Eq. (|B3|) . This assumption 
is indeed satisfied in our discussions on all configurations 
except ripples in tubules. Thus V 2 </> should be replaced 
with V 2 0- V • S in Eq. (JB3J) as well as Eq. (JTUl) when this 
condition is not met. 



2. Variation respect to the deformation of the 
surface 

In this subsection, let us denote Fq as the functional 
([§]) with h and kf vanishing, and define the additional 



F a d = 



k f 2 . 



dA. 



(B4) 



Any small deformation of a CLM without free edges 
can always be achieved from small normal displacement 
at each point r in the surface. That is, Sr = fl 3 e 3 . 
The frame is also changed because of the deformation of 
the surface, which is denoted as 



Se i = Q ij e j , (i=l,2,3), 



(B5) 



where fL 



-Qji, = 1,2,3) corresponds to the ro- 



tation of the frame due to the deformation of the surface. 
SF was first dealt with in Ref. [T3] which gives 

SF = J (2k c V 2 H - 2XH + P)fl 3 dA 

+ J k c {2H + c )(2H 2 - c H - 2^)^3^.(66) 

Fo llowin g Ref. [H|, and considering Eqs. (TATO]) , $M7\i 
and (|A18[) , through somewhat involved calculations, we 
can derive 

5{r m dA) = 2Qi 2 {H - n m )dA + V ■ mdfl 3 A m • dr 
+m ■ dr A d(dfl3 A *m ■ dr/dA), (B7) 

5(v 2 dA) = 2dQ 12 A *{d(f) + cv 12 ) 

+2dfl 3 A v • de 3 + 2{k v - H)v 2 fl 3 dA. (B8) 

Combining with the above two equations and using the 
integral by parts and Stokes' theorem, we derive 

SF ad = J h[V ■ (mV x m) + V x (mV • m)]n 3 dA 
+ J kj[{n v - H)w 2 - Vv: \7e 3 ]fl 3 dA, (B9) 

Then the Euler-Lagrange equation (|TT|) follows from 
SF = 5Fo + SF a d — because fl 3 is an arbitrary function. 



APPENDIX C: DERIVATION OF THE 
EULER-LAGRANGE EQUATIONS AND 
BOUNDARY CONDITIONS OF CLMS WITH 
FREE EDGES 

A CLM with an free edge can be described as a surface 
with an boundary curve shown in Fig. [TJ3. 



1. Variation respect to (f> 

Let 6(f) = S, we still have Eq. (|B3|) whose second term 
is transformed into 



kf <j) Vb^,ds 



(CI) 



with Vb — v ■ b. Thus from SF — we can derive the 
Euler-Lagrange equation (|33|) and boundary condition 
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Variation respect to the deformation of the 
surface 



Because one can select an arbitrary frame 
{r; ei, e2, 63}, here we take it such that ei and e 2 
align with t and b in the boundary curve. Then any 
small deformation of a CLM with an free edge can 
always be expressed as the linear superposition of small 
normal displacement ^3 and tangent displacement fl 2 
along e 2 at each point r in the surface. 

First, we consider the in- plane deformation mode 
6r = il 2 e 2 . The change of the frame is still denoted 
as Eq. (|B5|) . In terms of Rcf . .§£} , we have 5 § ds = 
— § Kg^ds where n g is the geodesic curvature of the 
boundary curve. Additionally, we can derive 



5 / GdA 



Gn 2 ds. 



(C2) 



Although the derivation of the above equation is some- 
what involved, its physical meaning is quite clear. Under 
the displacement f2 2 along e 2 , G is similar to the surface 
tension and the area element near the boundary decreases 
by {l 2 ds. Thus the variation of J GdA gives Eq. (|C2[) and 
then we obtain the boundary condition ([55)1 from SF = 0. 

Second, we consider the out-of-plane deformation 
mode 5r = £l 3 e 3 . Let us denote Fo as the functional 
(|32[) with h and kf vanishing, and define the additional 
functional as Eq. (|B4[) . SFq is fully discussed in Ref. [35[ 
as 

SF = J [k c (2H + c )(2iJ 2 - cqH - 2K) - 2\H}Q 3 dA 
2k c (V 2 H)VL 3 dA - f [k c (2H + c ) - kK n ]Q, 23 ds 



(C3) 



-2k c dH/dh + jn n — kT g ]£l 3 ds. 



Here the only difference is that the sign of k adopted in 
the present paper is opposite to that in Ref. [35| . Similar 
to the derivation of Eq. ([B9f from Eqs. (|B7|) and (fB8j) by 
using the integral by parts and Stokes' theorem, we have 

5F ad = J h[V • (mV x m) + V x (mV • m)]n 3 dA 

fe/[(« v - H)v 2 ~ Vv: Ve 3 ]tt 3 dA 

+ <p[h(vt + $) sin 2(f) — kfK n v t ]fl 3 ds 

■(h/2)sm24>n 23 ds, (C4) 



where Vt = v • t. In the above equations (|C3|) and (|C4[) . 
£l 3 represents the arbitrary small displacement of point 
on the surface along e3 and Q23 is the arbitrary small ro- 
tation of e3 around t at the edge. Thus SF = will give 
the Euler-Lagrange equation (|34[) and boundary condi- 
tions (|3"T)) and 



APPENDIX D: DERIVATION OF EQUATIONS 

flUD AND (fill) FOR TUBULES WITH 
NONUNIFORM TIL TING STATE IN SECTION 

[mD] 

For a cylindrical surface, take a frame such that ei, 
e 2 and e3 are along the circumferential, axial and radial 
directions, respectively. Then we have 



u)\ = ds, L02 = dz 
a = —l/p,b = c = Q. 



(Dl) 
(D2) 



Thus 



2H = a + c= -l/p,K = ac-b 2 = 0. (D3) 



Because dds — and ddz — 0, using the structure equa- 
tion (|A3|) . we have wi 2 = 0. Considering Eq. (|A14|1 . we 
obtain the spin connection 



S = 0, 

and then 

v = = 4> s ei + 4> z e 2 . 
From Eqs. (|A6|) and (|A13|) . we have 

K m = - COS 2 <j)/p, 



'88 ^ H'ZZ- 



(D4) 



(D5) 



(D6) 
(D7) 



Substituting the above two equations and Eq. (|E)3|) into 
the Euler-Lagrange equation (fTU)) , we can derive Eq. (TTJJ 
in Sec. ImDl 

Additionally, using Eqs. ^2]) . (|A4|) and (|D5|) . we de- 
rive 



v ■ de 3 = (cj) s /p)ds. 
From Eq. (|A19[) we can derive 

Vv: Ve 3 = -(pax/ p. 



(D8) 



(D9) 



In this derivation, one should note that dA — ds A dz = 
-dz A ds. From Eqs. flX20j) , ([MI]) . (jD2)) and ([D5]) . we 
can obtain 



(«v - ff)v 2 = (0 2 - 2 )/2p, 



(D10) 



and 



V • (mV x m) + V x (mV • m) (Dll) 

= 2sin2(/.(0 2 - <j) 2 s + cj) sz ) + cos2cj)((t) ss - <f) zz + A<j) z cf) s ). 



Substituting Eqs. (|D3| and (|D9|) - (|DTT|) into the Euler- 
Lagrange equation (llip . we can obtain Eq. (fT5)l in 
Sec. liLTTJl 
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APPENDIX E: DERIVATION OF EQUATIONS 
(pOjl AND (HU) FOR HE LICAL RIPPLES IN 
SECTION HUEl 

For the surface expressed as vector form 
r = Ml + y) cos(s/p),p(1 + y) sin(s/p), z} (El) 

with |y| <C 1, the first and second fundamental forms of 
the surface can be calculated as 

I = (1 + 2y)ds 2 + dz 2 , (E2) 

II = --(1 + y- p 2 y ss )ds 2 + 2py sz dsdz + py zz dz 2 (E3) 

P 

up to the order of O(y), respectively. The same order 
is kept in the following expressions in this section. In 
terms of the correspondence relations I = cof + u 2 and 
II = au> 2 + 2bujiuj2 + k>§> we have 

ui = (1 + y)ds, uj 2 = dz, (E4) 
a = -(l-y- p 2 y ss )/p, b = py sz , c = /n/ zz , (E5) 
2if = -(1 - y - p 2 y ss - P 2 y zz )/p, K = -y zz .(E6) 

2H is a function of s and z, using V 2 {2H)dA = d*d(2H), 
we can derive 



V 2 {2H) = (y, 



SS ~T UzZ 



P Ussss 



2p 2 y zzss + p 2 Vzzzz)/p 
(E7) 

Because du>i — y z dz Ads and du>2 = ddz = 0, using the 
structure equation (|A3|) . we have ui\2 = —y z oj\. Consid- 
ering Eqs. (|A14[) and (|A15|) . we have 

S=y z e 1 ,V-S = y zs . (E8) 

Using Eqs. (|aT2|) and (jATj]) . we obtain 

V0 = ^(l-»)ei + ^,ea (E9) 
V 2 = SS + 2Z + z y z - </> s y s - 2y0 ss . (E10) 



and then 

v = V0-S = [^,(l-|/)-p,]e 1 



(Ell) 



From Eqs. (|A2]l . ([ATf . (|Al9|) . ([Alp) , and (|ElT|) , we can 

derive 

Vv: Ve 3 = -0„/p + (y z <f) s + 2y^ sz + y zz )/p, (E12) 

(«v - H)V 2 = (0 2 _ 02 )/2 p + s (y0, + y 2 )/p 

+2py sz <t> s <t>z + (y + p 2 y ss - P 2 y Z z){(t> 2 s - <^)/2p.(Ei3) 

Since we only consider the case of h <C 1, up to the order 
0(/i), we obtain 

h(K m - H) = -(ft/2p)cos20, (E14) 
ft[V • (mV x m) + V x (mV • m)] (E15) 
= /j[2(0 2 - 2 + cj) sz ) sin 20 + (0 SS - 4> zz + 40 Z S ) cos 20] 



by using Eqs. (|A6]l and (|D11|) . 

Substituting the above equations (|E6[) - (|E15|) into the 
Euler-Lagrange equations (JTUJ) and (fTTj) . we can derive 
the basic equations (|20p and (l2~l]) under the assumption 
c = 0, A = 0, k f ~ fc c in SecUmi 



APPENDIX F: DERIVATION OF EQUA TIONS 
p9j) AND ([30]) FOR TORI IN SECTION IHTfI 



A torus can be expressed as a vector form 
r = {(r + pcosip) cos 9, (r + p cos ip) sin 6, psincp} (Fl) 
with r > p. A frame {ei, e2, 63} is taken as 
ei = {— sin#, cos 6*, 0}, 

e 2 = { — sin ip cos 9, — sin if sin 9, cos 99} , (F2) 
e.3 = {cos 9 cos sin# cos ip, sin 99}. 

By using Eqs. (|A1|) . (IF1[) and (|F2[1 . we can derive 

u>i = (r + p cos if)d9,u>2 = pdcp. (F3) 

From Eq. (|A2]l and (|F2|). we obtain 

W12 = [sin f/(r + p cos ip)]ui , (F4) 
W13 = -[cos<p/(r + pcos(/3)]wi, (F5) 

^23 = -(l/p)u> 2 . (F6) 

Considering Eqs. (|A14j) . (|A15|) and ijFi]). we have 

S = — {sin ip /(r + pcosip)]ei, (F7) 
V • S = 0. (F8) 



Additionally, (|A4[) . (TF5J) and ((F6J) give 

a = —cost/?/ (r + p cos f),b — 0, c = — 1/p. (F9) 



Thus by considering Eqs. ©, l(A5]) . (|A6]l . and (jAT3|) , we 
obtain 



?* + 2pcos(/5 cosy 



p(r + p cos if) p(r + pcosip) ' 



cos 2 4> cos <f sin 2 1 



r + p cos ip p 

r sin 2(f) 
2p(r + p cos if) ' 



(r + pcosip) 2 p 2 p(r + pcosip) 



(F10) 
(Fll) 
(F12) 
(F13) 



Substituting the above four equations into Eq. (|10[) . we 
can derive Eq. (J29J) in Sec. IIIIFI 
If = -7r/4, then V0 = 0, and 

v = \/(p— S = [sin if /(r + p cos (p)]ei, (F14) 
v • de^ = fiW3i = [sin2y/2(r + p cos if )]d9 ,(F 15) 

where v\ is the component of v in the direction of ei. 
Thus from Eq. (|A19[) we can derive 

Vv: Ve3 = — (r cos 2ip + pcos 3 ip)/p(r + pcos ip) 3 . (F16) 

Using Eqs. ([J20)) . ([Mi"]) . (|F9|) and (|FT4)1 . we can obtain 

(k v - H)v 2 = r sin 2 ip/ 2 p(r + pcosip) 3 , (F17) 
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and 



V • (mV x m) + V x (mV • m) = — cos <p/p(r + pcos<p). 

(F18) 

Additionally, from V 2 (2H)dA = d * d(2H), we have 



V 2 (2H) = r(r cos (p + p)/p 2 (r + pcosip) 3 



(F19) 



Substituting the above equations (|F10p and (|F16[> 
§FT9\i into Eq. dTTJ) , we can derive Eq. ([30]) in Sec. ITIlFl 



APPENDIX G: DERIVATION OF EQUATIONS 
(glD-(Il6]), AND flUJ) FOR T WIST ED RIBBONS IN 
SECTION HVCl 

A twisted ribbon can be expressed as a vector form 

r = {u cos ip, u sin 95, a(p} (Gl) 
A frame {ei,e2,e3} is taken as 
ei = {cos ip, sinip, 0}, 

e 2 = {— u simp, u cos ip,a}/\/u 2 + a 2 , (G2) 



e 3 = { a sin ip, —a cos ip, u} / \/u 2 + a 2 . 
By using Eqs. (|AT|l . (|G1~|) and ([G2j) . we can derive 



uji = du, u>2 = V u 2 + a 2 dip. 
From Eq. (fA"2j) and (|G2)) . we obtain 



(G3) 



W12 = kt/(u 2 + a 2 )]w 2 , (G4) 
wis = -[a/(u 2 + a 2 )]uj 2 , (G5) 
w 23 = -[a/(u 2 + e* 2 )H. (G6) 



Considering Eqs. (|A14j) . (|A15|) and ([G4]) . we have 



S = -[ U /(ii 2 + a 2 )]e 2! 
V • s = 0. 



(G7) 
(G8) 

Additionally, (|A4]l . ([G5]) and JG6| give 

a = c = 0,6=-a/(u 2 + a 2 ). (G9) 



Thus by considering Eqs. ©, fM)) . (|A"6)) , and (|A13|) . we 
obtain 



H = Q,K = -a 2 /{u 2 + a 2 ) 2 , 

n m — —asm2<p/(u 2 + a 2 ), 
.2,2 



(G10) 
(Gil) 

r m = -a cos 2^/(<it 2 + a 2 ), (G12) 
V 2 = 0„ tl + (w<?! )tl + w )/(u 2 +a 2 ). (G13) 



Substituting the above four equations into Eq. (|33[) . we 
can derive Eq. (J44J in Sec. IIVCI 
If 4> = or vr/2, then V0 = 0, and 

v = V0-S= [u/( U 2 + a 2 )]e 2 , (G14) 

v • de 3 = V2UJ32 = [ua/(u 2 + a 2 ) 2 ]du, (G15) 
where V2 is the component of v in the direction of e 2 . 
Thus from Eq. (|A19|) we can derive 



Vv: Ve 3 = 0. 



(G16) 



Considering sin 20 = and cos 20 = ±1, from 
Eqs. (|A2Q}, (jMH), ilSJ and (fCTI)) . we can obtain 

(rtv - i?)v 2 = 0, (G17) 
V • (mV x m) + V x (mV • m) = 0. (G18) 



Substituting the above equations (|G10P and (|G16j) ~ 
(jGl8l) into Eq. (O, we can derive Eq. gnj) in Sec. HVCl 

Now let us turn to the boundary conditions. If (f> = 
for a < 0, we have 4> = if - At the boundary, t||e 2 and 
u = u = W/2, so 



Vt = u / (uq + a 2 ), v b = 0. 
From Eqs. (|GI0)) - ((GT2|1 . we can obtain 



(G19) 



K = ~a 2 /(u 2 + a 2 ) 2 , (G20) 
K m = 0,T m = -a/{u 2 + a 2 ). (G21) 

Additionally, we have 

K n = acos 2 (7r/2) + &sin7r + c 2 sin(7r/2) = 0, (G22) 
Kg = ^12/^2 = uq/ (u 2 + a 2 ), (G23) 
T g = bcosn + (c - a)(sin7r)/2 = a/(ul + a 2 ),(G24) 

in terms of the definition of normal curvature, geodesic 
curvature and geodesic torsion. From the above equa- 
tions (|G19|) — (IG24|1 . we find that among the boundary 
conditions (|35|) (f38j) . only Eq. (|36| is nontrivial, which 
gives Eq. (|46|1 . Similarly, we obtain the same equation as 
(|46]) when </> = ir/2 for a > 0. 

The average energy per area (|48|) can be obtained from 
F/A, where A is the area of the twisted ribbon and F is 
the free energy (|32|) . The calculation is straightforward 
with c = 0, Eqs. fG10| . (ICT2]l . (iGTil) . and 



c?j4 = \J u 2 + a 2 dudip, 



ds = d Uq + a 2 dip. 



(G25) 
(G26) 
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